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Torques on interfaces can be described by a divergence-free tensor which is fully encoded in the 
geometry. This tensor consists of two terms, one originating in the couple of the stress, the other 
capturing an intrinsic contribution due to curvature. In analogy to the description of forces in 
terms of a stress tensor, the torque on a particle can be expressed as a line integral along a contour 
surrounding the particle. Interactions between particles mediated by a fluid membrane are studied 
within this framework. In particular, torque balance places a strong constraint on the shape of the 
membrane. Symmetric two-particle configurations admit simple analytical expressions which are 
valid in the fully nonlinear regime; in particular, the problem may be solved exactly in the case of 
two membrane-bound parallel cylinders. This apparently simple system provides some flavor of the 
remarkably subtle nonlinear behavior associated with membrane-mediated interactions. 

PACS numbers: 87.16.Dg, 68.03. Cd, 02.40.Hw 
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I. INTRODUCTION 

The crucial role of fluid membranes as a component 
of living cells was recognized early on by biologists 
In recent years, the mechanics of membranes has be- 
come a fertile area of research among physicists as well 
3|. Many cellular tasks such as exo- or endocytosis 
3|, the formation of vesicles 0, @, or the interac- 
tion with the cytoskeleton [tJ rely heavily on mechani- 
cal membrane properties and can be studied using the 
mathematical toolbox of theoretical physics. Various 
key questions are of an essentially geometrical nature: 
what is the shape adopted by a membrane subjected 
to some specified boundary conditions 8]? How are 
forces and torques transmitted through the membrane 
3, llOl. UA j . Il2l ? Which interactions does this imply 

m M MMM M M M ID, [11,1111? 

To answer questions such as these, one approach is to 
apply an effective formulation in which the membrane is 
modeled as a two-dimensional surface in space. This ap- 
proach is valid provided the length scales of interest are 
much longer than the membrane thickness which typi- 
cally is a few nanometers. The energy of the membrane 
is then completely determined by the geometry of the 
surface. In distinction to a simple fluid-fluid interface 
where the energy is dominated by surface tension [25j . 
bending elasticity now becomes important [III 1J- It 
is possible, in principle, to determine equilibrium pro- 
files by solving a field equation (the "shape equation") 
which follows from the minimization of the energy with 
respect to surface deformations. The forces and torques 
transmitted by the membrane are then imprinted on its 
shape. 

The shape equation is a fourth order nonlinear par- 
tial differential equation; as such, it can only be solved 



analytically in a few exceptional cases. One approach 
to sidestepping this obstacle is to expand the energy up 
to quadratic order. Unfortunately, this approach has the 
drawback that it is only valid if the surface is a weak per- 
turbation of some underlying reference shape, such as a 
plane or a sphere. If the membrane is highly curved, lin- 
earization becomes inadequate and other strategies need 
to be developed. 

In Ref. [j| an alternative approach was presented. It 
was shown that it is possible to relate the forces and 
torques to the local membrane geometry without any 
need to solve the shape equation itself. The link is formed 
by the surface stress tensor and its torque counterpart. 
This approach has already proven its value when applied 
to the problem of surface-mediated interactions [HI, [3, 
ll : particles bound to a membrane interact because they 
deform its shape [M M, EE 03 11 10, 11 11, 11 . In the 
traditional approach to this problem based on energy, the 
force between particles is calculated in three steps: first, 
the equilibrium shape of the membrane is determined 
for a given placement of the particles by solving the lin- 
earized shape equation. The energy of the surface is then 
evaluated by integration, and finally differentiated with 
respect to appropriate placement variables of the bound 
particles. This approach relies heavily on the validity of 
the linearization. Instead, one can also relate the force 
between the particles directly to the local geometry of the 
surface via the surface stress tensor. It becomes possi- 
ble to establish a host of valuable exact nonlinear results 

In this context the particles were fixed by horizontal 
forces allowing all other degrees of freedom, among them 
their vertical positions or tilts, to equilibrate. To gain 
further insight a generalization is necessary which will 
be provided here: we will extend the framework devel- 
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oped in Refs. [l3| and [lij and explicitly address ques- 
tions concerning the balance of torques. This will impose 
strong additional constraints on the shape of the mem- 
brane. To be concrete, we will examine symmetric two- 
particle configurations, thereby obtaining a transparent 
analytical expression describing the balance of torques 
which is valid in the fully nonlinear regime. We will ex- 
amine exactly the interaction of two membrane-bound 
parallel cylinders. This is possible because the membrane 
jrofile behaves like a one-dimensional elasticum of Eulcr 
28| , a wormlike chain in the language of polymer physics 
23 . l30l . l3l| ; thus it can be integrated exactly. This ap- 
parently simple system displays remarkably subtle be- 
havior. For strong membrane deformations the results 
of the small gradient approximation [2lj | fail and nonlin- 
ear effects, such as multiple solutions for one given set of 
parameters, emerge. Moreover, the simple exponential 
decay of the force or the cylinder tilt angle as a function 
of separation, which is correct for small deformations [2l| . 
gives way to a more complicated dependency. 

The paper is organized as follows: in Sec. [IT] we will dis- 
cuss how torques on surface patches can be determined 
within the framework of Refs. To begin with, we 

will treat a geometrical Hamiltonian which is given by 
a surface integral over a local density. The behavior of 
this Hamiltonian under boundary rotations will be used 
to identify the torque tensor. This tensor will be writ- 
ten down explicitly for a fluid bilayer. Surface- mediated 
interactions between particles on an asymptotically flat 
membrane are considered in Sec. IIII1 Total force and 
torque balance will be exploited to identify constraints 
on the geometry of the membrane. In Sec. IIVI the mem- 
brane geometry corresponding to two parallel cylinders 
will be determined exactly and the balance of forces and 
torques examined in detail. Finally, we summarize our 
main conclusions in Sec. |Vj 



II. TORQUES ON SURFACE PATCHES 
A. Surface energetics 

Consider an interface which can be modeled as a 
two-dimensional surface in three-dimensional Euclidean 
space. Such a surface E is described locally by its posi- 
tion X(^ 1 ,^ 2 ) € E 3 , where the £ a are any suitable set of 
local coordinates on the surface. The tangent vectors of 
E, e Q = dX/d£ a = d a X (a, b € {1, 2}), form a local co- 
ordinate frame; together with the extrinsic unit normal 
vector, n = e\ x eij\e\ x e^\, they define two geometri- 
cal tensors on the surface: the metric Qab — e a ■ and 
the extrinsic curvature K a b = e a ■ dbfi [32[ ■ The trace of 
the extrinsic curvature, K = g a bK ab , is the sum of the 
two principal curvatures whereas the Gaussian curvature 
Kq = det(i^) is their product. The covariant derivative 
on the surface will be denoted by V Q , the corresponding 
Laplacian by A = V Q V a . 

We associate with the surface an energy which can be 



described as a surface integral over a local density TL 

H S [X] = f dATL(g ab ,K ab ) , (1) 

where dA = ^det {g a b) d 2 £ = ^fgd 2 ^ is the infinitesi- 
mal area element. The density TL depends exclusively 
on surface scalars constructed using the metric and the 
extrinsic curvature tensor. 

To determine the equilibrium shape of the surface the 
functional ([1]) is minimized with respect to deformations 
of E. These deformations are described by a change in the 
embedding functions X — > X + SX. The corresponding 
variation of the Hamiltonian is [1(3] 

SH S = f dA E n-5X+ [ dA V a [H ab e b -Sn~f a - SX] . 

(2) 

Its bulk part is a surface integral over the Euler-Lagrange 
derivative £ (TL) times the normal projection of the sur- 
face variation SX. The second term is a surface integral 
over a divergence and can thus be recast as a boundary 
integral using the divergence theorem [33|]. It involves 
two parts: a part proportional to SX identifies the sur- 
face stress tensor f a [9|, [l(| ; the second part proportional 
to Sn involves TL ab = dTL/dK a b and originates from the 
partial integration of second derivatives of the deforma- 
tion. Note that additional terms would appear in the 
boundary term if TL also depended on derivatives of g ab 
and K ab . 

According to Noether's theorem, every continuous 
symmetry of the Hamiltonian implies a conserved cur- 
rent. The one associated with the translational invari- 
ance of H is the surface stress tensor f a . Its conserva- 
tion describes the balance of forces across curves drawn 
on the surface [1, [l(|. In an analogous way, the rota- 
tional invariance of H implies the existence of another 
conserved tensor m a which, as will be demonstrated in 
the next section, one can identify as the torque tensor. 



B. Identification of the torque tensor 

Consider a region of surface bounded by an outer limit- 
ing curve <9E out and suppose that the bulk of this surface 
is in equilibrium with N particles that are attached to it 
[3]. We now choose a curve around each of the particles, 
labeling these curves dT^i,i = 1, . . . , N. The outer curve 
<9E out together with all <9Ei encloses the surface patch 
E (see Fig. [1]). Considering this patch first allows us to 
use variation © to identify the torque tensor as no par- 
ticles are attached to E and every point on it is free to 
equilibrate. 

Under a constant infinitesimal rotation Set of one of the 
boundaries <9Ei, the position and normal vectors change 
according to SX = Sot x X and Sn = Sot x n. In this 
case, the boundary variation of the Hamiltonian (TTJ) is 
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FIG. 1: Surface region £ with 3 disjoint boundary compo- 
nents <9£i and an outer limiting boundary <9E out . 



equal to 

SH ® § _ §a . / dfi ^ [xx f a + H ab e b x n] 

= -6a- (f ds l^m a = -5a • M« xt , (3) 

JdY,i 



where the surface integral was converted into a boundary 
integral on the first line. The vector Z s = l^e a is the 
unit vector which is normal to <9£i and points out of 
the surface £; by construction it is tangential to X (see 
Fig. [J). The variable s measures the arclength along <9£i. 
In Eqn. ©, the boundary integral is identified as the 

external torque iVfjj ext acting on £ via its boundary dT^: 
the corresponding change in energy is given by (minus) 
the scalar product of M^ ext with the infinitesimal rota- 
tion angle 5a. The external torque M^ t on the surface 
patch £j is then simply given by — M^ ext due to torque 
balance. Notice that this argument does not require us 
to assume that the Euler-Lagrange equation is satisfied 
on Ej. 

To simplify notation we will restrict our considera- 
tions to patch without loss of generality and write 
M oxt := M^ t , with 



(4) 



M oxt = f ds l a m a , 

where I = l a e a = — i s . The surface tensor 
m a = X x f a + H ab e b x n 



(5) 



is the covariantly conserved torque tensor. It consists of 
two parts: a contribution due to the couple of the stress 
tensor f a about the origin as well as an intrinsic contri- 
bution proportional to Ti ab which originates only from 
curvature terms. The value of the former does depend 
on the choice of origin; however, its intrinsic counterpart 
is independent of this choice. One can easily check that 
the divergence of m a , V Q m a , is indeed zero in mechani- 
cal equilibrium as required by consistency with Noether's 
theorem. 



The occurence of curvature terms in the Hamiltonian 
is ultimately due to the fact that the physical surface is 
not infinitely thin but has an inhomogeneous force distri- 
bution along its transverse direction . If one captures 
this internal structure in the strictly two-dimensional 
Hamiltonian (|TJ) , the curvature terms give rise to an in- 
trinsic torque density and a non-vanishing normal com- 
ponent of the stress tensor. 

Note that the torque tensor, like the stress tensor, de- 
pends only on geometric properties of the surface (mod- 
ulo a contribution that is due to a shift of the origin). 
It is thus always possible in principle to determine the 
torques operating on a patch of surface from a knowl- 
edge of the shape of the surface alone. The snag is to 
determine the correct shape. 

It turns out, however, that exact geometric relations 
such as Eqns. ([4} and |5J are still very useful, even when 
we lack the details of the shape, if additional informa- 
tion, e^q. a symmetry, is available. This was shown in 
Refs. [HG3 in the context of surface-mediated interac- 
tions between particles where the focus of interest was 
the force. To see how this works for torques, let us first 
write down Eqn. ([5]) for the surface models we are inter- 
ested in: soap films and fluid membranes. 



C. Special Cases 

1. Soap film 

The Hamiltonian density of a soap film is given by 
Ti = er, where a is the constant surface tension. This 
implies that Ti ab — 0. Thus the intrinsic torque vanishes 
and, using f a = —oe a , one obtains @, 

m a = -a(X x e a ) . (6) 

The torque on a surface patch Si is therefore given by 



M cxt = -af ds (X x I) 



(7) 



It is simple to interpret this expression: the tangential 
stress — al provides a torque per unit length at every 
point of the contour <9£i. The line integral along <9£i 
yields the total external torque on the surface patch. 



2. Fluid lipid membrane 

A symmetric fluid membrane can be described by the 
Hamiltonian [26|, |27( 



H 



-K 2 + RKq , 



(8) 



where a is, as before, the surface tension, k the bending 
rigidity, and R the saddle-splay modulus. Together the 
two constants a and k provide a characteristic length 



A 



(9) 
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separating short length scales over which bending energy 
dominates from the large ones over which tension does. 

The last term of the Hamiltonian density (jSJ) can be 
written as the sum of a topological constant and a line 
integral over t he g eodesic curvature at the boundary of 
the membrane [321 ] . It does not contribute to the mem- 
brane stress tensor and does not enter the shape equa- 
tion. However, it gives a contribution to the torque tensor 
[34j : the derivative of the density ((H) with respect to K ab 
is 



H ab = nKg ab + R(Kg ab - K ab ) 



(10) 



Thus, the intrinsic torque does not vanish and, using the 
expression for f a from Ref. [3], the torque tensor can 
be written as 

m a = n{K ab - -Kg ab )K - ag ab (X x e b ) 
-K{V a K)(X x n) 

+ \(k + R)K g ab - RK ab \ (e b x n) . (11) 

Inserting this result into Eqn. yields the external 
torque 



£ ds[\^{Kl-Kj)-cj\{Xy.l) 

+ kK XI \K(X x t) - k(V±K)(X x n) 
t} , (12) 



nKt\ 



on the membrane patch Ei. In expression (|12p the unit 
tangent vector t = t a e a = n x I is introduced: it points 
along the integration contour i9Ei and is perpendicular to 
I and n. The projections of the extrinsic curvature onto 
the orthonormal basis of tangent vectors {I, i] are given 
by K± = l a l b K ab ,K n = t a t b K ab , and K ±[ \ = l a t b K ab . 
The symbol Vj_ = l a W a denotes the directional deriva- 
tive along the vector I. 

The intrinsic part of the torque tensor adds a contribu- 
tion to the integrand in Eqn. (| 1 2[) . which is proportional 
to the curvature K and tangential to the contour of in- 
tegration. Remarkably, the term M £ xt = J g ^ ds l a m% 
originating from the Gaussian curvature vanishes. The 
reason for this is that its integrand can be written as a 
derivative with respect to the arc-length s: 

l a m% = l a \R(Kg ab ~ K ab )] (e b x n) = -R(K\\t + K X{ \1) 



R(t a K ab e b ) = -itViin 



dn 

'"to' 



(13) 



where the Weingarten equations V a n = K e b were 
used and the symbol V|| = t a V a denotes the directional 
derivative along t, i. e. the derivative with respect to the 
arc-length. Integrating Vmti over the closed contour <9Ei 
yields zero and thus M £ xt = 0, even if an external torque 
is acting and even though the torque tensor comprises 
terms proportional to R. This result implies that torques 



on membrane patches do not depend on the last term of 
the Hamiltonian density ([5]) involving the Gaussian cur- 
vature Kq. As shape and stresses are also independent 
of that term, we can neglect it in the following consider- 
ations. 

In addition, if K is constant on <9Ei, the intrinsic part 
of the torque tensor does not contribute to M ex t at all. 
This is because the integral of t = dX /ds vanishes along 
any closed contour. 

It is straightforward to derive expressions analogous to 
Eqn. (JT3J) for any Hamiltonian of the form fl}. 

We will now show how the framework provided by the 
torque tensor may be applied to the study of membrane- 
mediated interactions. 



III. INTERFACE MEDIATED INTERACTIONS 
BETWEEN PARTICLES ON ASYMPTOTICALLY 
FLAT INTERFACES 

Particles bound to an interface may interact via the 
deformations they impose on the surface geometry. Con- 
sider, for instance, a multiparticle configuration such as 
the one depicted in Fig. [TJ forces as well as torques act 
on each of the particles due to the distortions in the 
shape caused by the others. Such a situation can thus 
only be stationary if the particle positions and orienta- 
tions are constrained by external forces and torques, re- 
spectively. As we have shown, these forces and torques 
are determined completely by the surface geometry. The 
external torque acting on the particle that adheres to 
the surface patch Ei, for instance, is equal to the exter- 
nal torque M cxt acting on the whole patch as given by 
Eqn. (|U); this is because the torque tensor is divergence- 
free on any part of the surface not acted upon externally. 
The torque we will consider in the following is the 

torque on the particle i mediated by the interface coun- 
teracting the corresponding external torque; the former 
is evidently equal to minus the latter (for i = 1 one thus 



has M 



(i) 



Similar arguments apply to the 
force F^' on the i th particle as explained in detail in 
Refs. IHIil. 



-M ext ) 
•th 



A. Balance of forces and torques 

The total force and the total torque must vanish if the 
multiparticle configuration is in mechanical equilibrium: 



N 



<i) - 



and 



JV 



M (i) = M c< 



(14) 



where -F ut is the force and iW ou t the torque on the outer 
boundary <9E out (see again Fig.[T]). 

In the following, we will consider, in particular, a sym- 
metric fluid membrane with nonvanishing surface ten- 
sion a characterized by the Hamiltonian density Q. Its 
ground state in the absence of particles is an infinite flat 
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FIG. 2: Two identical particles bound to an interface. As 
described in the text, it is possible to deform the contour of 
integration in order to exploit the available symmetries [13 ]. 



plane with zero curvature. In this case the energy is 
proportional to the area and thus infinite. This infinite 
constant in the energy will play no role in the sequel, and 
we will recalibrate the energy to set it to zero. 

Particles that are bound to the membrane typically 
deform its shape and increase its energy. These energy 
changes, however, must be finite whenever the forces and 
torques applied to fix the particle configuration are finite. 
This implies that the interface becomes asymptotically 
flat remote from the particles, even if a is infinitesimally 
small (see App. |A"|) . 

It is useful to decompose forces and torques into a hor- 
izontal part parallel to the asymptotic plane and a verti- 
cal part orthogonal to it. When the surface tension van- 
ishes, Kim et al. have shown that the vertical force and 
the horizontal torque on the outer boundary at infinity 
must vanish in equilibrium [23]. In App. E] we extend 
their discussion to treat situations in which a =/= 0. The 
change is non-trivial: it turns out that, whereas the ver- 
tical force vanishes as before, the horizontal torque does 
not necessarily vanish. In antisymmetric configurations, 
such as the ones that we will consider in the next section, 
it is this horizontal contribution on the outer boundary 
that makes it possible to balance an external torque. 

In the following section, the complete set of equations 
describing force and total torque balance [see Eqns. (fT4]) ] 
in a symmetrical two-particle configuration will be de- 
rived. 



B. Two- particle configurations with symmetry 

Consider a symmetrical configuration of two identi- 
cal particles bound to an asymptotically flat surface 
as described in Refs. [H H] (see Fig. W). We label 
by {x,y,z} the Cartesian orthonormal basis vectors of 
three-dimensional space K 3 adapted to the asymptotic 
plane. The vectors x and y lie parallel to this plane 
whereas z is its upward pointing unit normal. 

We will discuss two possible symmetries: mirror sym- 
metry in the (y, z) plane (the symmetric case) or a 
twofold symmetry with respect to the y axis (the an- 
tisymmetric case). The line joining corresponding points 
on the particles lies parallel to the (x, z) plane. We place 



the origin of the coordinate system in the middle between 
the two particles on the intersection line of asymptotic 
and symmetry plane (symmetric case) or the line of sym- 
metry (antisymmetric case), respectively. 

External torques on the surface can either be vertical 
or horizontal. The former cannot be compensated by 
the outer boundary. The physical explanation for this is 
that the surface does not resist shear forces. Horizontal 
torques, however, can be balanced (see App. which 
is necessary for a stable configuration: the separation 
between the particles is fixed by horizontal constraining 
forces. In the antisymmetric case, the two particles do 
not generally lie on the same line of action. Thus the 
forces will apply a horizontal torque M y to the surface 
which has to be balanced by the outer boundary. In the 
following, we will consider situations in which this torque 
is the only external torque on the entire surface. 

This restriction does not exclude external vertical 
torques on the individual particles; the symmetry 

will not be broken as long as all these torques cancel. 
Think, for instance, of a symmetric configuration consist- 
ing of two spheres on a soap film with a saddle-shaped 
(quadrupolar) line of contact (35| . A vertical torque M z 
on one of the particles is consistent with the symmetry so 
long as a torque —M z is applied to its partner. Such pos- 
sibilities may, in principle, be accommodated within the 
formalism. Here, however, we will only consider situa- 
tions in which the particle orientations have equilibrated 
and these torques vanish. 

As explained previously, the torque on the left particle 
can be obtained by integrating the appropriate projec- 
tion of the torque tensor along a contour surrounding 
the particle 

M = - f ds l a m a . (15) 

We choose the same contour of integration as used in 
Refs. [H|,[ll[ (see Fig. it consists of four branches with 
branch 1 lying on the symmetry plane/line and branches 
2-4 pulled open to infinity. 

Using Eqn. (| 14[) together with the assumptions made 
above we conclude that 

M = X p x F = My , (16) 

where F = Fx is the force on the left particle and X p 
is the position vector pointing from the fulcrum (i. e. the 
point about which the torque is acting) to the point where 
the force F is applied. We evaluate the torque about the 
origin in the Euclidean coordinate system {x, y, z}. The 
torque on the outer boundary enclosing the particles is 
then given by M out = 2M in the antisymmetric case; it 
vanishes in the symmetric case. 

The condition (Til)]) , together with the explicit expres- 
sion for the torque M derived here as well as its counter- 
part for the force F derived in Refs. [HI, [l4| , permit one 
to establish nonlinear relationships between geometrical 
quantities such as particle penetration or height differ- 
ences for a given model. 
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C. Equilibrium conditions from torque balance 



function 



1. 



General expressions for two finite-sized particles 



For a fluid membrane, the torque M on a particle can 
be obtained by inserting the corresponding torque ten- 
sor (fTTj) into Eqn. (jTHI) • As it evidently has to cancel the 
external torque on the particle, it can also be read off 
directly from the right hand side of Eqn. (p~2|) with the 
replacement of 9Si by the contour 1—4 and an additional 
minus sign out front. 

The resulting expression can be simplified further if 
one takes into account that the surface is asymptotically 
flat and possesses a certain symmetry. For finite-sized 
particles it is shown in App. [A] that the curvature terms 
of the torque tensor do not contribute to the torque on 
the membrane at infinity; the only term in the expres- 
sion for the torque due to bending stems from branch 1. 
Furthermore, the term proportional to (X x t) vanishes 
for the two symmetries because K±\\ K = on branch 1 
[3] . The torque M on the left particle is thus given by 



M = 



a f ds (X x I) — k / ds 

1-4 Ji 



- (V±K)(X xn)-Kt 
a j> ds (X x I) — k J ds 



\{Kl-Kl){Xxl) 



X x f - Kt . (17) 



where / = [\{K\ - Kf) - A~ 2 ] I - (V ± K) n. 

One now might naively expect that, in analogy to the 
curvature terms, the term due to surface tension vanishes 
along branches 2-4: the membrane becomes asymptoti- 
cally flat at infinity and the corresponding force has no 
vertical component. As X lies in the asymptotic plane 
due to our choice of origin, it would appear to be self- 
evident that the line integral vanishes. 

If this argument were correct, it would also apply to 
the complete outer boundary enclosing both particles and 
it would be impossible to balance the external torque. 
How then can a stable equilibrium be attained in the 
antisymmetric case if, as was argued earlier, an external 
torque is unavoidable? 

Fortunately, a closer inspection of the integral shows 
that it does not necessarily always vanish: for when 
branches 2-4 are sent to infinity, while the vertical force 
component indeed does converge to zero, the length \X\ 
of the corresponding lever arm goes to infinity simultane- 
ously. Additionally, one integrates along a contour whose 
length depends linearly on \X\. These two effects permit 
a finite value to remain: if we parametrize the remote 
regions of the membrane in terms of a height function 
above the asymptotic plane and consider its series ex- 
pansion (|A5j) for large \X\, a non- vanishing horizontal 
torque remains if the slowest decaying term is of order 
|-X"| _1 . In fact, from Eqn. (| A5|) we obtain for the height 



h(\X\,tp) = (Ci cos </? + £>i sin (^) |X| 
+ faster decaying terms 



(18) 



where C± and D± are constants and ip is the azimuthal 
angle measured from the x axis. Thus, branches 2-4 
can contribute a torque if either of the two constants in 
Eqn. (fT8|) does not vanish. A closer inspection of the 
line integral along branches 2-4 reveals that the term 
with coefficient C\ contributes to a torque about the y 
axis whereas the one including D\ does so about x (cf. 
calculation in App. [A] for M out ). The latter coefficient 
has to be equal to zero as we require the x component 
of the torque M out at the complete outer boundary to 
vanish (see again App. IA"|) . 

Bearing these considerations in mind, one obtains the 
following important relations by combining Eqns. (|16[) 
and (fT7|) . 



x-M = 
y-M = 

z-M = 



-k^cIs X-(fxx) = , 
X-(lxy) 

X-(fxy)-K(yt 
X-(fxz)-K(z-t) 



o~ (p ds 

J2-4 

—k J ds 



- k I ds 



(19a) 
(19b) 
F(z-X p ), 
= , (19c) 



where F = —k L ds (f ■ x) — a J 3 ds [l4|. These are the 
expressions referred to at the end of the previous sec- 
tion. They place strong constraints on the geometry of 
the membrane and can be simplified further for the two 
different symmetries. 

a. Symmetric case In the symmetric case, the 
derivative of K in the direction of I along branch 1, 
Vj_K, is zero. The curvature tensor is diagonal in (l,t) 
coordinates and thus vanishes. This implies that 

K\\t is equal to VyTZ [cf. Eqn. fig])], and ^dsK^t = 

^\y— oo ^|y— — oo 0* 

Furthermore, branches 2-4 do not provide a torque: 
the term in the height function (fTgj) that could give rise 
to a torque about the y axis is forbidden by the symmetry, 
i. e. its constant coefficient C\ = 0. Thus, 



a I ds (X x I) = 

J2-4 



(20) 



and Eqn (|17p simplifies to the expression 

M sym = - K Jds{ [^(Kl-K()-\- 2 ] (Xx^-iTj) , 

(21) 

which depends only on geometrical properties of the 
membrane at the symmetry plane. The force acting at 
each point along the mid-line is horizontal and Eqn. (|19ap 
is fulfilled identically. The two relations (|19b[) and (|19c|) 
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turn into 



« J^ds f sym (z ■ X) - Kj_(y ■ t) = F sym (z-X p ), 



(22) 



n J ds 



f Byui (y • X) + K x (z ■ t) 



= 



(23) 



with / s . 



[\{Kl - K 2 ) - A- 2 ] and F sym = oAL 



§ J 1 ds (Kj_ — KT) where AL is the excess length of 
branch 1 compared to branch 3 [T3 |. 

b. Antisymmetric case In the antisymmetric case 
branch 1 is a straight line and K± = Kii = 0. The 
torque is now given by 



M 



antisym 



g (p ds{Xxl) + K ds (ViK)(Ixn 



(24) 

In contrast to the symmetric case, the torque does not 
possess an internal part because the curvature K van- 
ishes everywhere along the contour. The force acting at 
each point along the mid-line is perpendicular to the y 
axis. As the origin lies on this line the only components 
of the torque transmitted through branch 1 which are po- 
tentially non- vanishing are the x and the z component. 
We have chosen, however, to set these components to 
zero [see previous Section and Eqns. (|19all9c[) ]. Thus, 
the contribution from branch 1 vanishes 

a J^ds {X x I) + n J^ds (V±K)(X x n) =0, (25) 

and the torque can be written as an integral along the 
boundaries 2-4 

cr / ds(Xxl) (26) 

J2-4 

where F antisym = k f t ds [\~ 2 (x ■ I - 1) + (V±K)(x ■ n)] 

This equation implies a strong constraint on the 
asymptotics of the membrane. Far from the particles 
the term of leading order in the height function (|18p has 
to be proportional to Otherwise, the external 

torque would not be balanced. The coefficient C\ has to 
be finite, its value fixed by the external force responsible 
for the torque. 

If the height function of the membrane is known, 
can be easily calculated: in App. [S] we have 



M 



antisym 



M 



antisym 



derived expression (|A7|) for M out in terms of the co- 
efficients of the height function (| A5|) / (TT5)) . From this 
expression we obtain Af ant i sym = (l/2)Af out = iraCiy, 
thus fixing C\ = M out / \2iro) . 



2. Two infinitely-long cylinders 

The situation changes if we consider two infinitely-long 
cylinders of radius R parallel to the y axis that are sep- 
arated by a distance d (see Fig. [3]). The size of the par- 
ticles is no longer finite. Thus, it is necessary to exercise 
a little care when carrying over expressions derived for 
finite-sized particles. 



a. Symmetric case The symmetric case corresponds 
to two cylinders adhering to one side of the membrane. 
The curvature parallel to branch 1, K\\, vanishes as 
branch 1 degenerates into a line. The contributions from 
branch 2 and 4 cancel. No torque is exerted from branch 
3: in contrast to the configurations with two finite-sized 
particles, the shape equation does not admit a term pro- 
portional to in the height function (see Sec. IIV|) . 

The torque per length L of the cylinder with respect to 
the origin of the Euclidean coordinate system {x, y, z} 
is thus given by 



M 



(27) 



where the length ho is the distance between the origin and 
the mid-point of the profile (see Fig. [3]). It is positive if 
this point lies below the origin. 

Let h c := — (z ■ X p ) be the vertical distance between 
the center of the cylinder and the asymptotic plane (see 
again Fig. [3]). It is positive if the two cylinders are located 
below the asymptotic plane. 

Torque balance establishes a relationship between ho 
and h c , involving only the geometry on the midline, 



-^sym,cyl 

where F s , 



given by Eqn ([27]): M 



-K 2 



sym,cyl 

Il4j implies 



c ^sym,cyl 



kK± - \K\h c 



K x - \K \h c 

x- 



±K 2 



(28) 



Eqn. (f2"5| remains valid if the interface is a soap film. 
The bending rigidity k then equals zero and we immedi- 
ately obtain the result that ho = 0. This is confirmed by 
the observation that the interface is flat everywhere and 
the forces on the left and right cylinder vanish [14j . 

If we parametrize the profile as a height h(x) above the 
asymptotic reference plane, the curvature K± at branch 
1 is exactly equal to —h"(0) in the symmetric geometry, 
where dashes denote derivatives with respect to x. In 
Ref. [2l| the height function h(x) was determined at the 
linearized level in the small gradient regime. Inserting it 
into Eqn. (j2"5|) validates this relation within the accuracy 
of the linear approximation (see App. IB)) . 

b. Antisymmetric case In the antisymmetric case 
the two cylinders adhere to opposite sides of the mem- 
brane. The torque on the left cylinder is determined as 
follows: at branch 1, the line integral is zero. The con- 
tributions from branches 2 and 4 cancel as before. Thus, 
the torque at branch 3 had better not vanish if the hori- 
zontal torque due the force F is to be balanced. However, 
the vertical force decreases faster than 1/x as x — > ±oo 
and thus cannot provide such a torque (c/. previous para- 
graph and Sec. lIV|) . How can torques balance under these 
circumstances? 

The solution to this apparent contradiction is the fol- 
lowing: the asymptotic plane for positive values of x does 
not necessarily coincide with the one for negative values. 
This is because the corresponding sections of the profile 
are disconnected. It becomes possible for the membrane 
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FIG. 3: Two parallel cylinders on a fluid membrane (symmetric and antisymmetric case). The profiles are determined as 
explained in Sec. II VI using the following parameters: a c — ai + a — 240°, d = 4, R = 1, and A = 1. The point about which 
the torques are evaluated is marked with o. 



to shift vertically with an offset /i ff at the origin with re- 
spect to the asymptotic plane(s) (see Fig. [3]). The torque 
is then simply given by 

-Mantisym,cyl/£ = ^ <7 ' 

This offset may be related to h c as follows. One 

has M ant i sym , C yi -(h c - % ffi )-F' a nti S ym,cyi where 

Fantisym.cyl/i = \J o 2 + (/cV ±K±) 2 - a Q- One obtains 



2h c ^1 + A 4 (V ± J^) 2 ' 

where h c :— [—(z ■ X p ) + h s/2] is now measured from 
the left asymptotic plane to the center of the cylinder 
(see again Fig. [3]). 

Expanding the inverse square root in Eqn. (|30[) up to 
zeroth and second order, respectively, yields a lower and 



an upper bound on the ratio of h s to h c : 

0<^< X 4 (y ± K ± f . (31) 

Obviously, the offset h g and h c have the same sign. If 
h c > 0, the left asymptotic plane lies above the right one 
as depicted in Fig. [31 for h c < the situation is reversed. 

By setting k — 0, we can again provide a check for 
consistency. The offset is zero as in the symmetric case. 
This is in agreement with the result from Ref. [l4[ that 
the interface is flat and no force acts on the cylinders. 

In Ref. [2l| the offset is set to zero in the ansatz for 
the height function. On first inspection, this may appear 
to be an error. However, it turns out to be consistent if 
Equation (|30[) is written in the small gradient regime: the 
first non-vanishing term is of second order in the small- 
ness parameter. As the height function of [2l| is itself 
correct only to first order, the two results agree at the 
level of the approximation. 

In this section several analytical conditions have been 
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derived which link different geometric properties of the 
interface profile to each other. For the antisymmetric 
case we have seen that the external torque is compen- 
sated either by a vertical force component (if the particles 
are finite) or by an offset (if the particles are infinitely- 
long cylinders). In the final part of this paper we will 
demonstrate the value of the framework developed in this 
section by showing how it can be applied to determine 
the exact shape of the membrane with two adhering cylin- 
ders. 



IV. TWO CYLINDERS ON A FLUID 
MEMBRANE - EXACT SOLUTION 

A. Determining the profile 

1. Shape equation 

If two parallel cylinders adhere to the membrane, the 
profile can be decomposed into the following parts: two 
bound sections in which the cylinder and membrane are 
in contact, an inner section between the cylinders, and 
two outer sections that become flat for x — > ±00 (see 
Fig. [3]). The shape of the bound parts is determined by 
the geometry of the attached particle, i. e. a circular arc; 
the profiles of the free membrane sections are determined 
by solving a nonlinear differential equation. Solving the 
equation itself is simple; the subtlety, as we will see, is in 
the implementation of the boundary conditions. 

We will introduce the "angle-arclength" parametriza- 
tion of the profile: the angle ip(s) between the x axis 
and the tangent to the profile as a function of arclength 
s completely describes the shape of the membrane (see 
Fig. [3]). It is connected to the curvature K via the re- 
lationship K = —ip, where the dot denotes a derivative 
with respect to s. 

The nonlinear shape equation that determines the pro- 
file of the free membrane is given by [1, [3(| 



2K + K 3 - 2\~ 2 K = 



It possesses the first integral 



K 



-K A - \- 2 K 2 = E 
4 



(32) 



(33) 



where E is a constant of integration. 

It is straightforward to integrate Eqn. (|3"3")l after a sep- 
aration of variables to determine s as a function of K 
which can be inverted to express if as a function of s. An 
additional integration yields the shape of the membrane, 
ip(s). One can do better however by taking advantage of 
the fact that the force -F me m per length L of the cylinder 
at every point of the membrane is not only constant but 
also horizontal on each membrane section. Thus, 



JL = -l a f a 



1 



-kK 2 



a ) I + kK n = /„ 



m x . 

(34) 



Projecting this equation onto n yields, with 77 = 

fmem/o- = COnst., 

X 2 ip — rysinT/i = . (35) 

where we have used x ■ n = — simp and K = —ip. 
Eqn. ([33)1 appears in a number of other physical applica- 
tions which include the motion of a simple pendulum 37] , 
the shape of a fluid meniscus under gravity [38j, or the 
behavior of elastic rods (Euler elastica) [28[ . Despite ap- 
pearances, it is completely equivalent to Eqn. (|33[) : both 
are of the same order in derivatives of ip; furthermore, as 
we will show, the constant of integration E can be written 
in terms of the scaled force rj and vice versa. Eqn. (|35p . 
however, has the advantage that it possesses a first in- 
tegral with an integration constant that can be deduced 
from Eqn. (f3"4"|) : the projection of the latter onto I yields, 
with x ■ I = cos ip, 



X 2 



-ip + 1] cos ip = 1 . 



(36) 



Eqn. (|3"6"j) can be interpreted as an energy balance of a 
fictitious particle moving in the potential V{ip) — rj cos?/), 
with displacement variable ip, mass A 2 and energy 1. The 
details of the "motion" depend on the value of ij (see 
below and Fig. |4j). From Eqn. (f36|) follows 



K = ±A"V 2 (l-?7CosV') • 



(37) 



It is simple to check that K satisfies the shape equa- 
tion (|32p . The constants n and E are not independent; 
the relation between them, r\ = zty/EX 4 + 1, is found by 
substituting K into Eqn. ([33]). 

In the following, we will determine the solutions of 
Eqn. (|3"6"j) for the free sections of the profile if the cylin- 
ders are separated by a distance d. Due to the symmetry 
it is sufficient to consider only the left half of the mem- 
brane (where x < 0, see Fig. While the outer sections 
are qualitatively identical for the two symmetries, the in- 
ner sections differ. Eqn. (|36[) is a first order differential 
equation involving one unknown constant r\. Thus two 
boundary conditions are required for each section. 



2. Outer section 

The boundary conditions for the outer section are the 
following: the free profile leaves the cylinder at a fixed 
contact angle ip(0) — a . We assume without loss of 
generality that a > and ip < [39| . At infinity, both 
ip and ip vanish. If we insert these two conditions into 
Eqn. p6p . we obtain rj — 1 on this section. This implies 
that E = everywhere along the profile. The corre- 
sponding potential V(ip) is plotted in Fig. 2] a): the so- 
lution with ip(0) = a Q decreases monotonically to ip = 
(where V = 1) as s — > 00. 

The corresponding shape equation can be solved in 
terms of elementary functions. One obtains 



i/j out (S) = 4arctan 



tan • 



(38) 
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FIG. 4: Above: Fictitious particle moving in the potential V(tp) = rj cost/' a) outer section b) inner section - symmetric case 
c) inner section - antisymmetric case; below: corresponding solutions (|38I40I42[1 of Eqn. (|36|l . 



where s :— s/\. Incidentally, the function also de- 
scribes the shape of a straight fluid meniscus that ap- 
proaches a solid surface at the angle a [1^, [38[ . 

The function ([38]) decays exponentially with increasing 
arclength. Thus, in contrast to the case of finite-sized 
particles, the horizontal torque due to a vertical force 
component, \M\ w cr\x\ smip out w cr|:z;|V> out vanishes at 
infinity (see App. [£]and Sec. IHIC|) . 



3. Inner section 

We first establish the connection between the total arc- 
length of the inner profile 2s m jd and the contact angle 
between the cylinder and the inner membrane, a\. If 
arc-length is measured from the mid-line, we have 

— — i?sino:i = / ds cos-0(s) , (39) 
2 Jo 

where all lengths, as before, are scaled with A. This con- 
dition determines the value of s m id implicitly in terms of 
a { . 

We again consider only positive values of ip. The so- 
lution for negative angles are obtained by reversing the 
sign of the functions (|4T)|) and (|4"2")) which corresponds to 
a reflection of the profile in the (x, y) plane. Note also 
that K = +ip for all shapes of the inner section due to 
the chosen orientation of the arclength. 



a. Symmetric case In the symmetric case, "0(0) is 
equal to zero. From Eqn. (|36|) we thus obtain rj < 1. 
For rj = 1, the cylinders are infinitely far apart and do 
not interact. If we omit this trivial case, the follo wing 
solutions for different values of ry are obtained t 28i [3l|, l40f 
(see Fig. H]b): 



7r + 2 am 

2am(W^ 
4 arctan 
arccos 1 



(gV^-K(m),m) , < 77 < 1 
, V = o 

1 < T] < 



, m 
tanh (| 



?7 < -1 

where am(s,m) is the Jacobi amplitude with parameter 
m and sn(s,m) = sin [am(s, m)]. The symbol K(m) de- 
notes the complete elliptic integral of the first kind [4l[ . 
The parameter m is given by m :— -^jL G [0, 2[. 

b. Antisymmetric case If the two cylinders adhere 
antisymmetrically, ip must vanish at the mid-line. Thus, 



cos ipj, 



(41) 



where "0mid = ip(0) is the angle at the mid-line. The case 
where ipmid equals zero or 180° corresponds again to the 
trivial solution with the two cylinders infinitely far apart. 
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No solution exists for 7/> mi d — 90° as the scaled force 77 
has to remain finite in an equilibrium situation. If < 
^ mid < 90°, 77 > 1 and, for -^K(i) < S < J=K(i) 
[ISIIH (see Fig. He) 



and 



antisym 



(~s) 



/ 2 

IS s SI 

I TO 



s^-K(I),l 



(42a) 

if it is greater than 90° and lower than 180°, rj < — 1 and, 
for - i K(i-) < a < -4=K(i)|l (see again Fig. 1 



antisym 



hi 

C 2 r 1 1 

(s) = arccos <1 sn 2 s-\/\r]\ — K( — ), — 

I to L m to J 

(42b) 

The parameter to is defined as above and varies now be- 
tween 1 and 2. 



4-. Boundary conditions at the cylinders 

To obtain the profile of the membrane for given sep- 
aration d and cylinder radius R, one has to determine 
the values of the scaled force 77 and the contact angles 
a Q and a\. The value of 77 for any given cti is determined 
implicitly by the requirement that 



^'"(smid) = a-i 



(43) 



where s m id is itself implicitly given by Condition (|39[) : the 
values of cti and a depend on the boundary conditions 
at the cylinder. We will discuss two possibilities here: we 
either fix the area of contact between the cylinders and 
the membrane, or consider attachment due to a finite 
adhesion energy 10 per area between the membrane and 
the cylinders. 

a. Fixed area of contact Suppose that the area of 
contact is fixed to 



a; = const 



(44) 



Torque balance will fix the relative values of ct and aj. 

The torque about the cylinder axis has to vanish in 
equilibrium. This is the case if the total energy E to t(cti) 
of the system exhibits a local extremum. The correspond- 
ing torque balance can be written as 



= Ki — K o — R(r) cos cti — cos a ) 



(45) 



where Ki := XKi and K a := XK Q . Eqn. (|45|) has two 
solutions: Ki = K Q and Ki — —K — 2/R. The latter 
implies that either Ki or K is smaller than — 1 / R. Since 
the cylinder has a curvature of — 1/R in scaled units and 
membrane and cylinder must not intersect, the second 
solution can be ruled out. The two contact curvatures 
must agree: 



dV> out 

Ki = K = 

as 



2sin(a /2) , 



(46) 



r\ cos a.\ = cos a 



(47) 



The values of 77, cti, eta, and s m id for given separation and 
cylinder radius can be determined numerically by solving 
the conditions (J39]), (g5J), (J44j) , and (gBJ) simultaneously. 
As an example, the two profiles for d = 4, R = 1, and 
a c = 240° are plotted in Fig. [3] . 

b. Adhesion balance If the cylinders attach to the 
membrane due to a finite adhesion energy w, the contact 
curvature condition holds (see [H, Sec. 12, problem 6], 
[H, and 0): 



Ki = K Q 



-1/R- 



(48) 



where w := . The curvatures are equal as in case 
a and torque balance (|45| is fulfilled automatically. The 
contact angle a can be determined from Eqn. (|46|) . Fi- 
nally, 77, Q!i, and s m id can be calculated using Eqns. (|39|) . 
(03J, and (jMl). 

The adhesion energy w and the wrapping angle a c are 
conjugate variables. Setting one of them to a constant 
value implies that the other will adjust in equilibrium. 
Questions of stability depend on which of the two vari- 
ables is fixed: a profile found to be stable under constant 
a c is not necessarily stable under the constant w. To 
avoid problems of this kind we will focus, in the follow- 
ing, on the constant wrapping angle scenario. 

Note that the boundary condition (|4"6")l does not always 
fix a unique profile. The total energy E io t{a.i) of the sys- 
tem may possess more than one minimum; thus profiles 
which are locally stable may exist. In the next paragraph 
we will see that this is indeed the case for certain ranges 
of values of d, R, and a c . 



B. Conditions from torque balance 

Now that we possess the complete profile we can go 
back and reexamine Eqns. (|28|) and ([30]) and study the 
behavior of ho := ho/X and h g := h g/X as a function 
of distance d for different values of R and a c . 

a. Symmetric case In the symmetric case, K± = 
y/2(l — 77) and Eqn. (|28|) can be written as 



where 



h = n- 1 [^2(l-7 1 )-{l-r ] )h c ] 



h c = 2sin(a Q /2) — i?cosa Q 



(49) 



(50) 



The height ho can be determined as a function of sepa- 
ration d for fixed cylinder radius R and different wrap- 
ping angles a c . If the value of a c is small, the solution 
coincides with the result of the small gradient approxi- 
mation (IB8j) as expected. The value of ho falls off ex- 
ponentially with the decay length I := 2 A for increasing 
separation d. 
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FIG. 5: The height ho at the mid-line as a function of separa- 
tion d for R — 1 and a c = 170° (long dashes), 180° (bold solid 
line), and 190° (short dashes). At d — 2 the two cylinders are 
in contact with each other and cannot draw nearer. For d < 3 
two locally stable solutions can be found which correspond to 
different angles a^. In the insets a and b the two solutions for 
d — 2.4 and a c = 170° are plotted where ai = 50.9° (a) and 
122.2° (b), respectively. Profile a corresponds to the global 
energy minimum here. 




-3.5- 



4.5 



► d 



FIG. 6: The tilt angle at (symmetric case) as a function of 
separation d for R = 1 and a c = 10° (short dashes), 45° 
(dashed-dotted line) , 90° (long dashes) , 180° (bold solid line) , 
and 240° (dotted line). 



For higher angles a c , however, the behavior changes 
at small values of d due to the breakdown of the small 
gradient approximation (see Fig. [5] for R = 1): the en- 
ergy Etot(cki) exhibits two minima instead of one [45[. 
These minima correspond to distinct stable profiles with 
different «i and ho- If a c = 180°, they possess the same 
energy. If a c < 180°, the global minimum is located 
at the smaller value of a; and the corresponding ho is 
always positive, i. e. the mid-point lies below the refer- 
ence plane (see again Fig. [5]). For larger contact angles 
a c this behavior is reversed: the higher value of a\ and 
the lower value of ho correspond to the global energy 
minimum. The mid-point may now even lie above the 
reference plane. 

In Fig. [S] the tilt angle at ■= %■ — cti is plotted as a 
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FIG. 7: The offset h g as a function of separation d for R — 1 
and different contact angles. The line styles are chosen as in 

Fig. ED 



function of distance d for R = 1. It decays exponentially 
with length A at large separations. If a c < 180° the tilt 
angle of the global energy minimum is always positive 
which corresponds to profiles such as the one plotted in 
Fig. a. If a c > 180°, the tilt angle of the global energy 
minimum is negative. For a c — 180° and for all d > 
3, at is equal to zero; at d w 3 the function bifurcates 
into two branches of same energy. A closer inspection 
of Eqn. (|47[) explains these findings: either a\ is exactly 
90° and thus at = 0, or r\ = — 1. At the bifurcation 
point both conditions hold. For d < 3 the profiles with 
r/ = — 1 are local minima of same energy whereas at = 
corresponds to the intermediate maximum we do not 
discuss further here. 

b. Antisymmetric case From Eqns. (|30p and (|35p 
one obtains with X7 ± K± = — i/>|s=o, K± = 0, and t] > 1 



Ks _ 1 _ 1 
2h c ~ V 



(51) 



where h c is given by Eqn. ([50)) . 

If the two cylinders are closer than 2i?, one additionally 
has to check for every stable profile whether the particles 
overlap or not. To avoid complications of this kind we 
will only consider separations d > 2R. In contrast to the 
symmetric case, one then finds a single solution for given 
d, i?, and a c . 

In Fig. [7] h g is plotted as a function of distance d. 
For small contact angles, h g is negative which implies 
that the asymptotic plane on the left lies below the one 
on the right and below the center of the left cylinder. 
For intermediate values of a c (such as 90° if R = 1) 
h s can be either positive or negative depending on the 
separation of the cylinders (see again Fig. [7]). If a c is 
further increased, h g is positive for all separations and 
the profiles resemble the one which is plotted in the lower 
part of Fig. [3] 

The tilt angle a t is always equal to zero for a c — 180°. 
The solution 77 = —1 which also solves Eqn. (TJ7J is now 
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FIG. 9: Scaled force r\ — 1 as a function of separation d for 
R = 1 and different contact angles (antisymmetric case) . The 
line styles are chosen as in Fig. H3 

forbidden due to the symmetry. For a c > 180° the tilt 
angle is positive; for a c < 180° it is negative. In both 
cases it decays exponentially with length A. 

C. Forces between the cylinders 

Using the stress tensor [IH, [lij], the force on the left 
cylinder is given by 

F cy i/L = f mem - a = a(rj- 1) . (52) 

Thus, one has only to determine the value of 77; the force 
follows directly. 

In Figs. [8] and [9] the absolute value of the scaled force, 
1 1] — 1|, is plotted for the symmetric and antisymmet- 
ric case. At large separations it decays exponentially as 
predicted by the linearized theory [21J. The underlying 
small gradient approximation, however, breaks down if 
the curvatures along the profile become too high, i. e. 
for large contact area and small separation. In Fig. [TU] 
the exact solutions are compared with the small gradient 
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FIG. 10: \r\ — 1| as a function of a c for d = 3 and R = 1 in the 
symmetric case. The small gradient approximation is plotted 
with a dashed line, the numerically exact solution with a solid 
line. 



approximation for the symmetric geometry by plotting 
1 77 — 1| as a function of a c for d = 3 and R= 1. The ap- 
proximate solution increasingly underestimates the force 
as the contact angle is increased. 

As soon as the small gradient approximation breaks 
down, different features appear: in the antisymmetric 
case, 1 77 — 1| increases faster for decreasing separation 
than a small gradient approximation would anticipate. 
In the symmetric case, two solutions can be found (c/. 
previous paragraph). The global energy minimum typi- 
cally corresponds to the solution with the smaller abso- 
lute value of the force. Remarkably, the force is constant 
[F cyl /(La) = -2] for a c = 180° and d < 3. 



V. CONCLUSIONS 

We have extended the formalism which was originally 
developed in Refs. 0, [l(| and primarily applied to forces 
[l3L [14] ] to study torques on surface patches. The un- 
derlying torque tensor depends only on geometric prop- 
erties of the surface (modulo the contribution that is 
due to a shift of the origin). Therefore, it is possible 
to relate torques on particles which are bound to an in- 
terface directly to the interface geometry. For the case 
of the Helfrich Hamiltonian we have explicitly worked 
out the torque tensor. Remarkably, its contribution from 
the Gaussian curvature part, even though nonzero, van- 
ishes upon closed-loop integration and therefore does not 
transmit torques between particles. 

By exploiting torque and force balance in the context of 
membrane-mediated interactions, several analytical con- 
ditions on the interface geometry have been derived for 
two-particle configurations obeying a certain symmetry. 
In the symmetric two-particle case all external torques on 
the whole surface can be set to zero; in contrast to that, 
such an external torque is inevitably applied in the anti- 
symmetric case due to the forces which fix the separation 
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between the particles. Total torque balance determines 
the asymptotic behavior of the membrane for both con- 
figurations. Particularly in the antisymmetric case, the 
torque is either compensated by a vertical force compo- 
nent (if the particles are finite) or by an offset between 
the two asymptotic planes (if the particles are infinitely- 
long cylinders) . Complications stemming from additional 
external torques ( e. g. due to a magnetic field if the two 
particles are dipoles) have not been considered here ex- 
plicitly but can be included with ease into the relevant 
expressions. 

To see how the rather abstract expressions from the 
first part can be applied, the two-cylinder case ("ID 
problem" ) has subsequently been solved in the nonlinear 
regime. The expressions obtained from torque balance 
could be used to calculate the values of the offset (anti- 
symmetric case) and the height at the mid-plane (sym- 
metric case). Finally, the force on the cylinders could 
be identified within the formalism as being essentially 
one of the parameters that had to be determined in the 
calculation of the profile. 

The simple ID problem already displays remarkably 
subtle behavior. For large deformations the small gra- 
dient approximation fails which is due to the intrinsic 
nonlinearity of the problem. The sign of the force, how- 
ever, is predicted correctly by the approximate solution. 
For two spheres on the membrane (2D problem), this 
need not be the case: although linear calculations pre- 
dict a repulsion for the symmetric configuration [20| , the 
nonlinear expression for the force F sym does not rule out 
that the interaction between the particles is attractive 
for higher deformations [cf. expression for F syul below 
Eqn. (|23p ]. While it is difficult to determine the sign an- 
alytically, coarsed-grained computer simulations suggest 
that a two-body attraction emerges if the deformations 
are high enough (irij . 
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APPENDIX A: FORCES AND TORQUES ON 
THE OUTER BOUNDARY 

Consider a symmetric fluid membrane with a finite 
number of finite-sized particles bound to it. What is the 
asymptotic behavior of profile properties such as the cur- 
vature K and its derivatives far away from the particles? 

To answer this question let us adapt a line of argument 
originally introduced in [23| for the case of vanishing sur- 



face tension. The excess membrane energy E induced by 
the particles is finite. From this requirement it immedi- 
ately follows that K has to be square integrable and thus 
vanishes far away from the particles. This implies that 
the behavior of the membrane is determined completely 
by surface tension there. Thus, the membrane becomes 
asymptotically flat. This fact allows us to use linearized 
Monge gauge to describe its behavior remote from the 
particles: the membrane profile is described in terms of 
its height h(r, ip) above the asymptotic plane where r and 
tp are the cylindrical coordinates defined on that plane; 
we choose to place the origin in the middle of the surface 
region to which the particles adhere. For large r, the 
curvature K is then equal to the negative (base plane) 
Laplacian of h, K = — V /i. 

The energy can be written as E = E a + E K + 
E R , where E a = (a/2) J dtp dr r (V/i) 2 and E K = 
(k/2) Jdipdrr (V 2 ft.) 2 . The contribution due to the 
Gaussian curvature, E R , is finite as it can be written 
as a topological constant plus a line integral over the 
geodesic curvature at the membrane boundary [32| . The 
energies E a and E K are both positive and thus finite as 
well. The corresponding membrane shape equation, ex- 
pressed in terms of K, is given by V 2 -K" — \~ 2 K which 
has the general solution [20( 

K(r,<p) = a K (r/X) 

oo 

+ y^(a„ cos nip + b n sin n<p) K n (r/X) 

n=l 

oo 

+ ^ (c„ cos rup + d n sin rvp) I n (r j A) , (Al) 

n=l 

where /„ and K n are modified Bessel functions of the 
first and second kind, respectively [4l| . The functions /„ 
tend to infinity as r — > oo. Thus, c„ and d n must vanish. 
The asymptotic expansion of K n for large r is given by 

and shows that the functions K n (r) decay essentially ex- 
ponentially for increasing r. 

The force -F ut on the outer boundary is a line inte- 
gral over the appropriate projection of the stress tensor, 
ids t^-slaf 1 - Since l a f a depends on derivatives and 
powers of K plus a surface tension term [l4| , the contri- 
bution to the force due to curvature vanishes as r — > oo. 
The only surviving contribution at infinity is a horizontal 
pulling force proportional to a. 

Although the vertical force vanishes at infinity, an ex- 
ternal torque M on t may still act on the outer membrane 
boundary. The reason for this is that the length of the 
lever arm increases linearly with R. The origin of the 
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torque is the surface tension and it is given by 

M out = a / dtp R[Rr + h(R, ip)z] x I 
Jo 

= a dtp R{-h(R,ip)l v r + [h(R,tp)l r - Rl z ]f 
Jo 

+Rl v z} , (A3) 

where I = l r r + l v <p> + l z z in cylindrical coordinates with 
basis vectors {r, cp, z}. The unit vectors r and ip lie 
on the reference plane whereas z is normal to it. The 
moving trihedron of the boundary curve consists of the 
vectors /, f , and n which are perpendicular to each other. 
In particular, I may be expanded in terms of the height 
function h and its derivatives: 



I 



k[(9A)(drh)}\ r=R <p + ±(d r h)\ r=R ; 



qR 



\l 1 + {iwii \ r=R \ + {^d r h)\ r=R } 

(A4) 

where q := 1 + [ ( ^ |r = R l . 

The height function h(r, tp) satisfies the linearized 
shape equation V 2 (V 2 /i) = A~ 2 V 2 /i. The solution of 
this equation consistent with a finite value for E a is [2(| 

h(r,<p) = C +A K Q (r/X) 

oo 

+ ~y^^(A n cos rup + B n sin rup)K n (r/ A) 

n=l 
oo 

+ (C n cos rup + D n sin nip) r~ n . ( A5) 

n=l 

Inserting Eqn. (|A4[) and (|A5|) into the equation of the 
torque (|A3|) one obtains 

""7i cosy? + A sin / 1 xi 

s V+ °li^JJ ' 

(A6) 

The unit vector <p still depends on the coordinate ip. To 
evaluate the integral over ip in (|A6[) . we rewrite <p as 
tp = — sin a; + cos y> y, where x and y are the Cartesian 
coordinates parallel to the asymptotic plane. We finally 
obtain for R —> oo 



M out = a I d<pR 
lo 



M out = 2na{-D lX + dy) 



(A7) 



Hence, the horizontal torque is in general not equal to 
zero even though the vertical force is. Its contributions 
stem from those terms of the height function (|A5[) which 
are proportional to r _1 and decay slowest for increasing 
distance from the particles. Conversly, the existence of 
an external torque forces the surface to display a very 
slow 1/r decay. 



APPENDIX B: CALCULATIONS ON THE 
LINEARIZED LEVEL 



In this appendix we will show that Eqn. (|28p is consis- 
tent with the results of Ref. [21[ . To this end we calculate 
h(0) — —ho and compare it with the right hand side of 



Due to the symmetry it is sufficient to consider only 
the left cylinder. Its axis is located at x = —d/2 (see 
Fig. [3]). For -d/2 - S a < x < -d/2 + 5 h cylinder and 
membrane are in contact (Si j0 = i?sinai j0 ). In small 
gradient approximation the shape of the outer membrane 
segment (x < d/2 — S a ) is given by [2l| 



/W (a;) = Sexp (-\x + d/2\/X) , 



(Bl) 



where B = —S X/R to first order in S Q (strictly speaking, 
the expansion is in 5 /X here and below). The membrane 
segment adhering to the left cylinder has the circular 
profile 



h cy \(x) 



h c + ^i? 2 ~{x + d/2 f 



(B2) 



where h c is defined in Fig. [3] The profile has to be con- 
tinuous at x = — d/2 — S . This condition yields 



-Vr^ 

-R 
R 



51 



A 

R ! 



SoX 
R 

0(5 2 



exp (— <S /A) 



Acoth^- „ . „, 



(B3) 



In the last step we exploited the fact that 6 = 
coth(^) Si to first order (2lj . 

The profile of the membrane between the two cylinders 
in the symmetric case is given by pH 



h in (x) = C + D cosh (x/X) 



(B4) 



where C and D are constants that can be determined 
from the conditions of continuous profile and slope at 
x = — d/2 + S[. The latter condition yields 



D 



XSi 



1 



XSi 



^W^s 

the former 
C = -{h c 



sinh(^r^) 



i?sinh(^) 



O(Sf), 
(B5) 



R 2 — S 2 ) — D cosh 



= -h c -R 



Acoth(^) 



R 



S i + 0(S 2 ) ^ 0{5 2 ) . (B6) 



To first order the depth h n at the mid-line of the profile 
[47l ] is thus given by 



ho = -MO) = -(C+D) 



XSi 



i?sinh(^) 



It can also be obtained from Eqn. 

-KM = -j? y ields 



-0{8 2 ) . (B7) 
Using K± = 



ho = - - 



D 



2A 7 



XSi 



1 



D 2 
2A 7 



i?sinh(^) 



0(5f). (B8) 



The two results coincide at first order which confirms the 
validity of Eqn. (|28[) within this approximation. 
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